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Abstract 

Dark energy candidates for which the equation-of-state parameter w is less than —1 violate the 
dominant energy condition, and are typically unstable. In scalar-tensor theories of gravity, how- 
ever, the expansion of the universe can mimic the behavior of general relativity with u; < — 1 
dark energy, without violating any energy conditions. We examine whether this possibility is phe- 
nomenologically viable by studying Brans-Dicke models and characterizing both the naturalness of 
the models themselves, and additional observational constraints from limits on the time-dependence 
of Newton's constant. We find that only highly contrived models would lead observers to measure 
w < -1. 
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I. INTRODUCTION 



The last several years have seen a sustained flow of cosmological data, beginning with 
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Dservations of type la supernovae bolstered by large-scale redshift surveys 



and measurements of the cosmic microwave background (CMB) from the ground 



the o 

!: .noo. mm . .... ..... . WM^P 

satellite j^. These studies have indicated that the expansion of the universe is accelerating 
(the scale factor obeys d > 0), and that the total amount of clustered matter in the universe 
is insufficient to account for the small value of its average spatial curvature. 

Both of these features (acceleration and flatness) can be explained in the context of 
convention al g eneral relativity by invoking a smooth, persistent dark-energy component. 

To be compatible with the observed isotropy and 
homogeneity of our universe on large scales, the energy-momentum tensor of the dark energy 
should be that of a perfect fluid. 



+ Pxgiiu , [i-j 

where is the fluid rest-frame four- velocity, px is the energy density and px is the pressure. 
The dark energy must be smoothly distributed in order to escape detection in the dynamics 
of gravitationally bound systems and large-scale structure. To make the universe accelerate, 
general relativity implies that the pressure px must be appreciable and negative. From the 
Friedmann equations 
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^H' = ^p-- (2) 



and 
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d • o AttG , . , , 

- = H + 2H' = — p + 3p , 3 
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we find that acceleration only occurs if the total pressure p (dominated by the dark energy, 
since matter is pressureless) is less than — p/3 (including both matter and dark energy). 

Since CMB observations imply that the spatial curvature is small (|K/a^| <^ H^), it is con- 
venient to set K = and consider a two-parameter set of cosmological models, characterized 
by the density parameter in matter = STTCpm/SH"^ and the dark energy equation-of-state 
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parameter wx, 



Through the continuity equation 



Wx = Px/px ■ 



px = -3H{1 + wx)px 



(4) 



(5) 



the equation-of-state parameter governs the rate at which the dark energy evolves as the 
universe expands, 

-3(1 + (6) 



din a 

A strictly constant vacuum energy (a cosmological constant) would have wx = — 1- Obser- 
vational constraints are often presented as exclusion contours in the VL^-wx plane. 

Of course, it is impossible in principle to directly measure the pressure in a component 
that is smoothly distributed, since there are no pressure gradients (or at least no readily 
observable ones). Constraints on wx actually derive from observations of the behavior of the 
scale factor a(t), and use the Friedmann equations ()2I3|) to translate these into limits on wx- 
Different observational methods will be sensitive to different integrated behaviors of the dark 
energy density; it is nevertheless useful to consider the instantaneous effective equation-of- 
state parameter that would be derived from ()2I3|) in a flat universe dominated by matter 
and dark energy: 

Weff ■ 



(7) 



where 

is the ratio of ener; 
bounds 



a 



VLx 1 — 

density in matter to that in dark energy. Current observational 
on yield 



1.48 < WcfT < -0.72 



(9) 



at the 95% confidence level. 

While the task of identifying a compelling candidate source with equation of state param- 
eter in this region is a formidable challenge for fundamental physics, the portion satisfying 
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Wes < — 1 is particularly troublesome theoretically. (Note that recent data, e.g. 
seem to suggest, at one or two sigma confidence, a Wes changing with z to values less than 
-1 today.) It is only possible to obtain w < — 1 by violating the dominant energy condition 
(DEC), which for a perfect fluid can be stated as 



P > \P\ 



(10) 



The physical motivation for the DEC is to prevent instability of the vacuum or propagation 



of energy outside the 



for some time 
by CaldweU 
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i ght cone. Nevertheless, such energy components have been known 



L gnt 
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ications o 



27| and their role as possible dark energy candidates was raised 
who referred to DEC-violating sources as "phantom" components. The 
;hese have since been investigated by several authors (for some examples see 
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Typically, DEC-violating 



sources with w < —1 are subject to violent instabilities, alth oug h these may conceivably be 
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cured in models with higher- derivative kinetic terms 

Despite the difficulties of model-building, it is certainly worthwhile to consider the 
possibility that Wes < — 1 when characterizing observational constraints, if only to keep 
open the possibility of a surprising discovery. In addition, given how little we under- 
stand about dark energy, we should keep an open mind about the true explanation for 
the acceleration of the universe. One alternative that has been explored is a modi- 
ication of gen eral relativity that would only become important on cosmological scales 

~^ 42,13 Q, y, 3,^2, y, y, 3, Q. riedmann equations ()2p3|l are not 



valid, the reconstruction (jZj) of the equation-of-state parameter would similarly be invalid. 
Therefore, it becomes conceivable that we could measure the effective value Wcs (as de- 
fined by (|7j) to be less than —1, even if the actual dark energy source obeys the DEC, or 
if there is indeed no dark energy at all. In other words, we could be tricked into thinking 
that Wx < — 1 by assuming the validity of general relativity (see also 
non-phantom matter behaving as phantom matter). 
/^From ((Zj) it is clear that we infer Wcs < — 1 today if 



for examples of 



where the subscript indicates present-day values of quantities. 

In this paper we consider a simple class of modifications of general relativity: scalar- 
tensor theories, featuring a scalar field that interacts non-minimally with gravitation {e.g., 
through a direct coupling to curvature). A wide variety of alternatives to GR can be cast 
as scalar-tensor theories, at least in some range of validity. Focusing on Brans-Dicke (BD) 
theories js^l, we examine whether a scalar-tensor theory of gravity could account for the 
acceleration of the universe and yield Wes < —1, while remaining consistent with other 
experimental constraints. In particular, the cosmological evolution of the scalar (p leads 
to time-dependence of Newton's gravitational constant G, which are constrained by solar- 
system tests of gravity. We find that only extremely unnatural and contrived models lead 
to an inference oi w < —1, even in this wide class of extensions to general relativity. 



II. BRANS-DICKE THEORIES 



A particular class of theories of gravity beyond General Relativity are Brans-Dicke theo- 
ries j57| with a potential, which may arise, for example, from the dimensional reduction of a 
higher dimensional theory. These theories consist of a metric and a Brans-Dicke scalar field 
(f, with action 



S- 



BD 



d X 
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ipR--{d^ip)d^ip-2V{ip) 



+ J (fx^f^CMii'i.g) , (12) 



where CM{'4'iy9) is the Lagrangian for matter fields ipi and t^; is a constant. In this frame, 
the Jordan frame, matter is minimally-coupled to gravity and hence test particles fall freely 
along geodesies of the metric g^i,. Scalar-tensor theories such as these have been considered 



previously as a way to solve the the fine-tuning problems of quintessence (see |59l l60l l6ll|.) 

The predictions of Brans-Dicke theories differ from those of GR due to the presence 
of a new scalar component to gravity. Since GR is well-tested in the solar system, these 
deviations must be smaller than the accuracy of current observations. This can happen in 
one of two ways. Either the potential is sufficiently confining to render the Brans-Dicke scalar 
essentially constant in the solar system or, in the presence of a sufficiently weak potential, 
the parameter u must satisfy the bound u > 40000, obtained using signal timing from the 



Cassini spacecraf t 16211 . However, this bound may be weaker on cosmological scales than in 
the solar system 63|. 

We will be interested in two possibilities. First, that a smoothly distributed background 
matter component, described by Cm and minimally coupled to the Brans-Dicke sector, may 
fuel faster than exponential expansion. Second, that the Brans-Dicke scalar itself may lead 
to a similar effect, without appealing to sources outside the gravitational sector of the theory. 

It is convenient to perform both a conformal transformation and a field redefinition to 
obtain an Einstein frame description of the theory. 

We define a canonically-normalized version of the Brans-Dicke scalar via 



where 



0"* = 



167rG0 , 



/2t<; + 3 



(13) 



(14) 



with uj the Brans-Dicke parameter. Denoting the Jordan- frame metric by g^^, the Einstein- 
frame metric is 

g^, = e^/'^'g,, (15) 

and the resulting action becomes 



St; 



—R-.(9,.a)rT.-Uia) 



+ / d^Xy^CM{'ipi,g,(T) 



(16) 



where 

U{a) = 2e-2'^/'^*F[v9(a)] . (17) 

This description now has the advantage that the gravitational sector is of Einstein form, with 
a minimally coupled scalar field, but the disadvantage that the matter does not freely-fall 
along the geodesies of the Einstein-frame metric. 

Further, the relationship between the energy-momentum tensors in the two frames is 

Therefore, if the energy conditions are not violated in the Einstein frame they will not be 
violated in the Jordan frame, since Q is strictly positive. 
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Let us now study the Einstein-frame equations of motion. We define the Einstein-frame 
scale factor and time coordinate by 

a = e"/2-*a 

dt = e"/2"*rft (19) 

and let a prime denote differentiation with respect to t, so that the Einstein-frame Hubble 
parameter is 

H=-. (20) 

ct 

The Friedmann equations are then 



— r, 871G , . , . 

H' = ^ (p. + p) , (21) 



H' = -A-rG 
and the equation of motion for a is 



a'Y + p + p , (22) 



a" + 3iy(T' + f/,, = -^(p-3p) . (23) 

Here we have also defined p = e^'^l^^p and p = 6^°"/°"*^, where p and p are respectively the 
Jordan frame energy density and pressure of matter. Since a is a canonically-normalized, 
minimally-coupled scalar field, its energy density and pressure are given by the usual defini- 
tions 

P. = + , (24) 
P. = - U{a) . (25) 

Because cosmological observations, and in particular those quantities entering (jllj) . in- 
volve Jordan frame quantities, we need to know how to transform between our easily- 
interpretable Einstein-frame quantities and those in the Jordan frame. This is achieved 
by the following expressions 
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and 



a" a' 



a 



H' - — + —\H 



(27) 



Substituting the Einstein-frame equations of motion into the expressions and 
we obtain, after some algebra 

SttG 



H 



{2UJ + 3) 

/ 47rG 



cj + 2) p + (cr 



J\2 



+ 



\2uj + 3 



and 



'\2 



2cj + 3 



(28) 



(29) 



These are the fundamental expressions relating the observable quantities H and H to sources 
as measured in the Einstein frame. 



III. A PERTURBATIVE APPROACH 

A. Acceleration from a Dark Energy Component in Brans-Dicke Theories 

Let us now focus on the case in which the source of cosmic acceleration is some dark-energy 
component in the matter sector (rather than the BD field itself). Since we are interested in 
the present day, we have a = ctq = 0. We assume that the BD field provides a negligible 
contribution to the evolution of the universe, so that 

(a')^ Uia) « p , (30) 

and that the theory is consistent with solar system test of gravity, requiring uj ^ 1. However, 
we will keep quantities that are first-order in \j ^/To^ since these are important for our effect. 
The "matter" sector will consist of ordinary matter plus the dark-energy component X. 
Expanding and evaluating (j^Hj) and ()29|) at the present day, after some algebra we obtain 

ijo^-47rG(l + ^-4e-r/)po , (31) 



8 



Hi ^ At^G ('^ - Po , (32) 



where we have introduced the parameters 
and 
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3(2cj + 3)po 



(33) 



V = I • (34) 

'167rG(2cj + 3)po 



(Note that, since = 0, we have po = Pxo + Pmo, so we could drop the tildes, but we won't 
bother.) We thus have 

I 2-3^ 

which implies 



^ -3 ( o/ ' 1 . (35) 



WeS ~ ^^7^^^ • (36) 

So far we have not assumed that the parameters ^ and r] are small. However, note that 

3(2cj + 3)po ■ ^ ' 

Since we have already assumed that (o"q)^ ^ po, and also that u; ^ 1, we see that ^ is 
certainly very small. It is therefore legitimate to rewrite (j36|l as 

Wes^wx- (l + ^<y]^-{l + a)r] . (38) 



Note that for a = 3/7, this is 



2 



29 10 

WeS ^ Wx - - ' (39) 



and since ^ is extremely small (less than 10~^), the second term on the right hand side will 
have a negligible effect on Wes- 

Now let us turn to rj. If we Taylor-expand the potential to first order about the present 
day value cr = 0, 



a 



U{a)^Uo{l + X-) , (40) 
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where f/(0) = Uq and U^^{0) = XUq/g^, we obtain 

2At/|, 

" " ■ 

Again, since we have assumed f/o ^ po and ^ 1, this looks very small except for the 
freedom associated with the dimensionless parameter A. Thus, it is possible to obtain 
detectably below —1 if 

X>ur^) . (42) 



vf/o/ 

Thus, we certainly can get Wcs < — 1 with matter sources that obey the Dominant Energy 
Condition. However, we should think about what it means to have rj be large while ^ remains 
small. From the definitions (jHHj) and (jS^), this occurs only if the time derivative of the BD 
field is small while the slope of its potential is very large. This will only happen for very 
finely-tuned conditions - either a very sudden change in the slope of the potential, or when 
the initial conditions for the field are chosen such that it has recently been climbing up 
the potential and is near a local maximum today. Hence, although it is possible to obtain 
Weff < — 1 behavior in BD theories, it is by no means natural. 



B. Acceleration from the Brans-Dicke Scalar with no Dark Energy 



The right hand side of the Einstein-frame Friedmann equations are sourced not only 
by matter, but also by the BD scalar a. It is natural to ask whether super-exponential 
acceleration can be obtained purely from this modification of general relativity. 

Making the approximation pUjl . using (jSEI), evaluated today and defining x = a'^j y/Pcr,o 
we obtain 

4 V2w + 3 



Ho 



3 2?i; + 4 
2 2^7 



■ X 



73 2w + 4 



^x^ + 



X + 



f/n 3A 



-3 



w + 2 



2Uo 

Pcrfl 

a 



2a 
l + a 



Pcr,o 2w + 3 2w + 3 1 + a 
The Friedmann equation can then be rewritten as 



Pcrfi 



1 + 



Aw + 6 



■ X 



X 

~2 



a 
l + a 



(43) 



(44) 
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which allows us to obtain 



Ho _ 1 [ 3 ^ 



—pz y/Aw + 6 — 2wx 



+ 3A-^-3-^(w + 2)| . (45) 

Pcr,o 1 + a I 



From equation (j44p we can see that in the limit -C 1, which we may term the slow-roll 
regime, we obtain 



' (46) 



which, from (^3)1 yields 



(y 

Wcff < -1 for A > - (47) 
iio > for A > a (a; + 2) 

Therefore, given that u is typically required to be much larger than unity, there exists a 
range in which the universe is not superaccelerating, but in which we infer w^s < even 
without an additional dark energy component. Again, though, we require a very large value 
of the parameter A characterizing the slope of the potential. 

IV. CONSTRAINTS FROM THE TIME- VARIATION OF NEWTON'S CON- 
STANT 

Thus far we have focused only on the behavior of the effective equation of state parameter. 
However, in Brans-Dicke theories it is important to remember that, for a weak BD potential, 
the field (f plays the role of a dynamical Newton constant G^s via 

Experimental constraints on the time variation of Newton's constant yield 

\GcS,o\ ^ in-12 
CcfT.O 

Since (jlSjl gives 



< 6 X 10"'^ yr"' . (49) 
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this then imphes 

< 6 X 10-^2 yr"^ , (51) 



or 



/2m; + 3 1 IGeffol /2w + 3 5.88 X 10-2 

^^^-^—ir.^^i— —' <^2> 

where we have written Hq = lOO/i km s^^Mpc"^, with h = 0.72 ± 0.08. 

We would hke to know what range of Wcs is possible today while remaining consistent 
with bounds on the time- variation of Newton's constant. To this end it is convenient to 
define 

y = jr r — • 

In terms of this variable equation becomes 

= liy) =1 + 7: 77^ y - ^ — ' (54) 



Pcr,o V 2 J 12 " 1 + a 

and using equation ()45|1 we obtain 



WcfT = -(1 + a) 



1 + o (4 - wy) + - — — 7(y) 



3'' ' 2w + 3 1 + a 2w + 3 



(55) 



Clearly, is a function of three parameters: Q^/^x = a = 3/7, A, giving the slope of 
the potential today, and the Brans-Dicke parameter u. For fixed u = 40000, we plot this 
expression as a function of the variables y and A in fig. ^ As can be seen, negative values of 
A correspond to positive Wes- 

Fixing the parameters A and u and plotting this function for \y\ < 5.88 x 10"^, we see 
that, in order to get values of Wcs ~ —1-1, we need A > 0.8 for = 1 (see fig.E)), or A > 5000 
for u = 40000 (see fig. El)- In this last case we must also demand \y\ < 0.001. 

From the graphs we see that in order to obtain Wcs ~ —1.1 with uj = 1 we require 
A ~ a. However, this condition implies a severe fine-tuning. The equation for the scalar 
field requires that the second derivative be negative and much larger in magnitude than the 
first derivative. Therefore, a must be close to a maximum today. This is consistent with the 
toy potential we used as a test. 
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50000 



1 

FIG. 1: u)eff as a function of y = -^f^ and A 



V. STUDY OF GENERAL SOLUTIONS FOR BD ACCELERATION 

To establish the main result of this paper, it was sufficient to use the perturbative approach 
of the previous two sections. However, it is interesting to consider the equations of motion 
(12 ()22|) and ()23|) in the Einstein frame in generality. Making the following definitions 



■7" — a/ Pcrfl t 

d = d/do 
U 



7 Pcrfi f — 



(56) 
(57) 
(58) 



and redefining a prime to denote a derivative with respect to r, the equations become 



SttG 



^ a 
2 1 + a 



(59) 



13 



-1.05 



-1.15 



-1.2 



-0.04 -0.02 



0.02 0.04 

y 



FIG. 2: Wef£ as a function of y = ^ with u = 1. 
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,'2 



a 



1 



+ a 2a^ 



q; 



(60) 
(61) 



We focus on the case of A > a/2 in order to ensure < —1 today. As we shall 
see, for some potentials (for example square or exponential ones), (t{t) has almost zero first 
derivative {x <^1) but a much larger and negative second derivative, impling that a achieves 
a local maximum today. This in turn implies a significant fine tuning of the initial conditions 
for a. 

One possibility is to impose that a' is always small, for example a' = x aX all times, so 
that a" ^ 0, and to look for the potential that would satisfy this condition. It is convenient 
to think of /(cr) = /[c(^)], so that / becomes a function of time. For example we could use 
the Priedmann equation, the scalar field equation and the second Einstein equation. Since 



1 da 



\/Pcr,o - -r 
^ a dr 



1 da 
d di 



Hq = Hq + 



^0 
2(7. 



(62) 
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FIG. 3: Wes as a function of y = ^ ^ with oj = 40000 and A » q. 
consistent initial conditions are 

a(0) = 1 (63) 

ct(0) = (65) 
/(O) = 1 . (66) 

We must then solve for the three new unknowns: a, cr, /. In this case we find that Wes is 
very close to —1 at all times, as shown in figure (j3)). 

An alternative possibility is to look for a potential / that leads to a constant less 
than -1. This requires 

£ - ^' 2^ - f'^ - - - T^: ■ 

Using this relation, the second Einstein equation and the scalar field equations, with 
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FIG. 4: WeS as a function of both 2/ = and A with lo = 40000. 
initial conditions 

a(0) = 1 

a'(0) = ^J^ + 
cr(0) = 
a'(0) = a; 
/(O) = 1 , 

we find that the potential is not even a function any more (/ becomes double- valued) . 
The potential is a curve and a = is a point at which the curve is continuous but not 
differentiable. This tells us something about how unnatural a potential like this is. So in 
this constant Wes does not seem likely at all. 

The general issue here is the following. Since cr' = x ^ 1, equation ()44|) implies that 
7 ~ 1/(1-1- a). Thus, in order to achieve < —1 (e.g. = —1-2) with u = 1, 
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(68) 

(69) 
(70) 
(71) 
(72) 



equation (jlTj) yields A ~ a. Further, equation (jH^ implies Hq ^ Hq/ y/Pcr,Q and therefore 
since a* i/o ~ 1; the scalar field equation (plj) evaluated at the present day is 

a 



° 2(1 + a) 
where we have used /,o-(0) = A/a*. This requires 



an 



> 1 or 



Mp > 1 . 



(73) 



(74) 



In other words, the function cr(r) must either be close to a maximum, having a small first 
derivative and a negative second derivative or, as in section IIIA, close to a sudden change 
in the potential. Similar behavior also occurs for uj = 40000 with A ~ 5000. 

VI. CONCLUSIONS 



The observed acceleration of the universe presents a fascinating yet daunting challenge 
to particle cosmology. The development of a theoretical framework in which the correct 
magnitude of the dark energy density and its relatively recent dominance are explained 
in a manner consistent with both particle physics and general relativity has so far eluded 
researchers. 

Another problem is posed by the equation of state parameter of the dark energy or, 
equivalently, the time evolution of its energy density. Even for the most economical possi- 
bility, a cosmo logical constant, we are faced with a problem that has confounded theorists 
for decades. The cosmological constant problem remains even when we turn to other dark 
energy models. U —1/3 > w > —1 then, assuming a solution to the cosmological constant 
problem, there are many proposals to obtain the requisite acceleration. To date, however, 
all of these face fine tuning problems or worse, when considered as serious particle physics 
or gravitational models. 

A more severe problem exists in the observationallv allowed ran ge w < — 1. Such a source 



xists m tne oDservationaiiv aiiowea ran ge w < — i. buci 

Q, y, 3, 0, H Q, H 0, Q, HQ y , 0, 



for the Einstein equations [la, |2a, 123, i2a 1221, laa laii laa, laa laa, laa |4^ |4l| must 

violate the dominant energy condition and hence may lead to instabilities of the vacuum or 
the propagation of energy outside the light cone. 
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One way around this problem is to consider theories in which the Friedmann equations 
are not vahd so that it is possible to measure the effective value Wes to be less than —1, even 
if the actual dark energy source obeys the DEC, or if there is no dark energy at all. 

In this paper we have considered Brans-Dicke (BD) theories, to see if such an effect 
can occur naturally in this framework, while remaining consistent with other experimental 
constraints such as those on the time- dependence of Newton's gravitational constant G. We 
have demonstrated that only extremely unnatural and contrived models lead to an inference 
of w < —1, even in this wide class of extensions to general relativity. In particular, it is 
necessary to fine tune the behavior of the Brans-Dicke scalar so that it is approaching a 
maximum today, having small first derivative with respect to time and yet large second 
derivative. 
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